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Abstract 


The Wide Field Survey Telescope (WFST) is a dedicated time-domain multi-band (u, g, r, i, and z) photometric 
survey facility under construction. In this paper, we present a preliminary study that assesses the quality of 
photometric redshifts based on WFST by utilizing mock observations derived with the galaxy catalog in the 
COSMOS/UltraVISTA field. We apply the template fitting technique to estimate photometric redshifts by using 
the ZEBRA photometric-redshift code and adopting a modified set of adaptive templates. We evaluate the bias 
(median relative offset between the output photometric redshifts and input redshifts), normalized median absolute 
deviation (Nmap) and outlier fraction (foutier) of photometric redshifts in two typical WFST observational cases, 
the single 30s exposure observations (hereafter shallow mode) and co-added 50 minutes exposure observations 
(hereafter deep mode). We find bias <0.006, onmap < 0.03, and foutier SS 5% in the shallow mode and bias 
0.005, onmap © 0.06, and foutiier ^: 1796-2796 in the deep mode, respectively, under various lunar phases. 
Combining the WFST mock observational data with that from the upcoming CSST and Euclid surveys, we 
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demonstrate that the Zphot results can be significantly improved, with foutlier ~ 1% and anman ~% 0.02. 
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1. Introduction 


The development of modern astronomy has given rise to an 
increasing demand for powerful multi-band photometric sky 
surveys. Such surveys, e.g., the Sloan Digital Sky Survey 
(SDSS; e.g., Brescia et al. 2014; Albareti et al. 2017; Zhao 
et al. 2021), Dark Energy Survey (DES; e.g., Drinkwater et al. 
2010; DES-Collaboration et al. 2016; Ivezié et al. 2019), and 
Hyper Suprime-Cam Subaru Strategic Program Survey (HSC- 
SSP; e.g., Aihara et al. 2018; Hikage et al. 2019), with well- 
designed equipments, reasonable observational strategies, and 
fruitful scientific results in stellar physics, galaxy physics, and 
cosmology, have demonstrated their strong impacts on modern 
astronomy. 

The Wide Field Survey Telescope (WFST) is a dedicated 
time-domain multi-band (u, g, r, i, and z) photometric survey 
facility under construction jointly by University of Science and 
Technology of China and Purple Mountain Observatory, which 
is expected to start commissioning observations around 2023 
August. WFST has a 2.5 m primary mirror, an active optical 
system, and a 0.73 Gigapixel mosaic CCD camera on the main 
focus plane; moreover, WFST is located near the summit of the 
Saishiteng Mountain in the Lenghu area that is a world-class 
observational site (Deng et al. 2021), thereby achieving high- 


5 Xiao-Zhi Lin and Yu Liu contributed equally to this work. 


quality imaging over a field of view of 6.5 deg”. The main 
science goals of WFST surveys are time-domain sciences 
including supernovae, tidal disruption events, multi-messenger 
events, and active galactic nuclei (AGNSs), asteroids and the 
solar system, the Milky Way and its satellite dwarf galaxies, 
and galaxy formation and cosmology (Wang et al. 2023). 
Robust determination of cosmological redshifts is one of the 
most crucial factors in fulfilling the above WFST science goals. 
However, high-precision galaxy redshift measurements require 
spectroscopic observations for each source (ie., obtaining 
spectroscopic redshifts, Zspec). This task is not only expensive 
but also time consuming. Alternatively, there is another way to 
measure redshifts using photometric surveys (i.e., obtaining 
photometric redshifts, Zphot), which is much more efficient than 
spectroscopic observations. This method, although not as 
precise as the z,,,, measurement, has demonstrated its 
extensive use in the Z,,.¢ determination of an overall huge 
amount of survey targets at one time (e.g., Benjamin et al. 
2010; Brescia et al. 2014; Cavuoti et al. 2017; Sánchez & 
Bernstein 2019). The application of Zphot has enabled a wide 
range of exciting extragalactic sciences as mentioned above. 
To date, a series of methods have been developed to estimate 
Zphot- In general, they can be divided into two main categories. 
One is based on template fitting that works as follows: the 
observed photometry is compared to a given set of pre-assumed 
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galaxy templates to determine the best-fit redshift corresp- 
onding to the maximum likelihood (e.g., Benítez 2000; 
Feldmann et al. 2006; Brammer et al. 2008; Luo et al. 2010; 
Rafferty et al. 2011; Yang et al. 2014; Cao et al. 2018). The 
other is the so-called training-set method, which constructs a 
neural network (e.g., Collister & Lahav 2004; Blake et al. 
2007; Sánchez et al. 2014; Pasquet et al. 2019) and performs 
machine learning to obtain Zphot, focusing on finding empirical 
relations between the redshift and galaxy properties (e.g., 
magnitudes and colors). This method is usually based on a 
large sample of secure Zspec, which are mostly available in the 
lower-redshift universe. However, since the magnitude limits of 
all WFST bands are deeper than most of the current Zspec 
surveys, it is difficult to find a sample of well-measured Zspec 
that can be representative of the full survey sample. Therefore, 
in this paper, we choose to measure Zphot of mock WFST 
observations based on the former technique, i.e., template 
fitting. 

The main goal of this paper is to assess Zpnot quality of the 
WFST photometry system preliminarily. We utilize the 
COSMOS /UltraVISTA multiwavelength galaxy photometry 
catalog (Muzzin et al. 2013) to produce mock WFST data. This 
survey has deeper magnitude limits than WFST observations, 
so it is suitable to select a subsample of galaxies from this 
survey, the magnitudes of which meet the WFST detection 
limits. Utilizing this subsample, we generate the mock flux of 
each WFST filter passband based on WFST instrumental 
parameters with good data quality, and then estimate the 
corresponding observational error. We choose to use the 
ZEBRA code (Feldmann et al. 2006) for Zpnot estimation. The 
main advantage of this code is that it can generate a new set of 
templates adaptive to the observations to minimize the 
mismatch between observed spectral energy distributions 
(SEDs) and the galaxy templates that are either from theoretical 
synthesis models or observed certain types of galaxy SEDs in 
the local universe, thereby improving Zphot quality. 

This paper is organized as follows. In Section 2, we 
introduce the WFST photometry system, COSMOS /Ultra- 
VISTA galaxy catalog, and generation of mock WFST data; in 
Section 3, we introduce the process of Zphot computation; in 
Section 4, we show Zpnot results and make comparisons with 
other works; and in Section 5, we summarize our results. All 
the magnitudes quoted are AB magnitudes. 


2. Data 
2.1. Overview of the WFST Photometry System 


WEST has six filters, i.e., u, g, r, i, z (see Figure 1) and w, 
with the white-light w band specifically designed for detecting 
asteroids in the solar system and thus being excluded from Zphot 
computation in this paper. There are two planned key programs 
of the 6 yr WEST survey: the wide-field survey (WFS) program 
and the deep high-cadence u-band survey (DHS) program. The 
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Figure 1. Total transmission curves of the WFST ugriz filters (indicated by 
different colors), with the instrument response and atmosphere absorption and 
scattering taken into account. 


WES program aims to survey a total of ~8000 deg? sky area in 
the u, g, r, and i bands in the northern hemisphere, with about 
90 visits in each band over 6 yr given a single exposure of 30s 
for each visit; while the DHS program plans to routinely 
monitor a total of 2:720 deg? sky area in the highly sensitive u 
band surrounding the equator every year, with a much higher 
observing cadence (down to hours) and being supplemented by 
a multi-band ancillary survey. The z-band imaging is excluded 
in the WFS program due to its relatively low efficiency and 
limited contribution to time-domain sciences; moreover, high- 
quality z-band imaging data will be achieved by other northern- 
hemisphere surveys such as Wide Imaging with Subaru HSC of 
the Euclid Sky (WISHES). However, WFST will allocate some 
additional observational time (about 1300hr over 6 yr) for 
specific purposes or particular interests, e.g., capturing time- 
critical targets and mapping the Galactic plane, which require 
intensive scanning of certain sky areas using the z-band 
imaging. 

In this paper, we compute Zphot in two typical WFST 
observational cases, i.e., the single 30 s exposure observations 
(hereafter shallow mode) and co-added 50 minutes exposure 
observations (hereafter deep mode). The deep mode can be 
realized by integrating all the observational time in each band 
mainly with the WFS program, thus achieving deeper detection 
limits than any existing single-telescope surveys with compar- 
able survey areas in the northern hemisphere (Lei et al. 2023; 
Wang et al. 2023). 

The average night sky background brightness at the WFST 
site (i.e., the Saishiteng Mountain, Lenghu Town, Qinghai 
Province) is approximately V = 22.0 mag arcsec ? when the 
moon is below the horizon; under new moon conditions, the 
best sky level can reach 22.3 mag arcsec ^? , which is measured 
in the extreme case when the bright part of the Galactic Disk is 
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Table 1 
5c Limiting Magnitudes and Sky Backgrounds of WFST Observations (from Lei et al. 2023) 
Lunar Phase Observational Mode WFST-u WFST-g WFST-r WFST-i WFST-z 
0 deg (no moon) Shallow mode 22.31 23.42 22.95 22.43 21.50 
Deep mode 24.86 25.95 25.48 24.96 24.03 
Sky background 23.27 22.82 21.80 20.99 20.05 
45 deg (1/4 moon) Shallow mode 22.27 23.30 22.89 22.40 21.49 
Deep mode 24.82 25.84 25.42 24.93 24.02 
Sky background 23.02 22.49 21.66 20.93 20.03 
90 deg (half moon) Shallow mode 22.04 22.86 22.62 22.26 21.43 
Deep mode 24.58 25.38 25.14 24.78 23.96 
Sky background 22.00 21.37 20.99 20.61 19.9 
135 deg (3/4 moon) Shallow mode 21.64 22.34 22.21 21.99 21.31 
Deep mode 24.17 24.85 24.72 24.51 23.83 
Sky background 20.86 20.21 20.08 20.01 19.61 
180 deg (full moon) Shallow mode 20.97 21.62 21.58 21.49 21.00 
Deep mode 23.48 24.12 24.09 24.01 23.51 
Sky background 19.30 18.73 18.78 18.97 18.92 


far away from the local zenith (Deng et al. 2021). In this paper, 
we set the sky background to a fixed value of 
V = 22.3 mag arcsec~*. Under this circumstance and with no 
moon, the 5c limiting magnitudes can reach depths of 
ugriz = [22.31, 23.42, 22.95, 22.43, 21.50] in the shallow 
mode and ugriz — [24.86, 25.95, 25.48, 24.96, 24.03] in the 
deep mode, respectively (Wang et al. 2023). The modeling 
results of the 5c limiting magnitudes and sky backgrounds in 
different lunar phases are listed in Table 1 (Lei et al. 2023). 


2.2. The COSMOS/UItraVISTA Galaxy Catalog 


In this paper, we adopt the multiwavelength galaxy 
photometry catalog in the COSMOS/UItraVISTA field 
(Muzzin et al. 2013) to produce mock WFST data, given 
that it has deep optical coverage, broadband photometry, and 
high-quality Zphot and corresponding best-fit galaxy SEDs. 

This catalog covers a sky area of 1.62 deg? with point-spread 
function (PSF) matched photometry in 30 bands, with the 
wavelength range extending from 0.15 to 24 um, including two 
ultraviolet bands (FUV and NUV) from the GALEX satellite 
(Martin et al. 2005), seven broadband (u, gu r'di yz, B; Vj) 
and 12 medium-band ([A427-IA827) optical data from the 
Subaru and Canada-France-Hawaii Telescope (Capak et al. 
2007; Taniguchi et al. 2007), four near-infrared imaging bands 
(Y, J, H, K;) from the UItraVISTA survey (McCracken et al. 
2012), and the 3.6, 4.5, 5.8, 8.0, and 24 um channels from 
Spitzer's IRAC and MIPS cameras (Sanders et al. 2007). The 
5c depths of the COSMOS/UltraVISTA survey in all bands 
are tabulated in Table 2, with typical depths in optical bands 
being deeper than those of WFST (see Table 1). 

Photometric redshifts of galaxies in the COSMOS /Ultra- 
VISTA catalog are computed based on the template-fitting 


technique with the EAZY photometric-redshift code (Brammer 
et al. 2008). The default seven EAZY templates are comprised 
of six templates derived from the PEGASE models (Fioc & 
Rocca-Volmerange 1999) and a red galaxy template from the 
models of Maraston (2005). To improve the quality of the 
fitting, Muzzin et al. (2013) added two additional galaxy 
templates: one is a one-gigayear-old single-burst galaxy 
template generated from the Bruzual & Charlot (2003) model 
to improve the template fitting for galaxies at z > 1 with post 
starburst-like features; and the other is a slightly dust-reddened 
young galaxy template to improve the fitting of UV bright 
Lyman break galaxies (LBGs) with heavy dust extinction at 
1.5 « z « 3.5. EAZY fits the observed multiwavelength photo- 
metry of galaxies utilizing linear combinations of the above 
nine initial templates (as shown in Figure 2) based on the x? 
minimization algorithm. Muzzin et al. (2013) provided in their 
COSMOS /UItraVISTA catalog the best template combination 
coefficients for each of the galaxies, so that we can generate its 
best-fit SED. We show some of the best-fit galaxy SED 
examples at their respective redshifts from the COSMOS/ 
UItraVISTA catalog in Figure 3. Photometric redshifts derived 
by Muzzin et al. (2013) are of high quality, being consistent 
With Zspec from the ZCOSMOS survey: up to z 1.5, their Zphot 
are accurate to Az/(1 + z) = 0.013, with an outlier fraction of 
only 1.6%; up to z~ 3, their Zphot show good agreement with 
Zphot from the NEWFIRM Medium Band Survey. 


2.3. Generation of Mock WFST Data 


First, the mock flux in each band for each galaxy in the given 
catalog can be calculated by convolving the galaxy redshifted 
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Figure 2. The set of nine initial galaxy templates adopted by both Muzzin et al. (2013) and this work for zs; derivation (the templates have been normalized 


appropriately for displaying purpose; see main texts for details). 


Table 2 

Depths of the 30 Bands in the COSMOS /UltraVISTA Photometry Catalog 
Filter FUV NUV u B; m V; rt it 
50 Depth 25.2 25.1 26.4 27.3 27.0 26.6 26.8 26.2 
Filter gr IA427 IA464 IA484 IA505 IA527 IA574 IA624 
So Depth 252 25.8 25.6 25.9 25.6 25.7 25.4 25.7 
Filter IA679 IA709 IA738 IA767 IA827 Y J H 
So Depth 25.3 25.4 25.4 25.1 25.1 24.6 24.4 23.9 
Filter K, 3.6 um 4.5 um 5.8 um 8.0 um 24 um 
So Depth 23.7 23.9 23.3 21.3 21.0 45 py 


SED with the filter transmission curve, using 
[2 $ARO)4A 


poss = = : 
[2 M0)2A 


(1) 


where S, is the best-fit observed SED of the COSMOS/ 
UltraVISTA galaxy, and R(A) is the transmission curve of one 
of the five WEST filters. The mock flux F?"*** is then calibrated 
to the mock observational flux according to the i band apparent 
magnitude (Subaru i* flux, a) given in the COSMOS/ 
UltraVISTA galaxy catalog. This conversion is performed by 
using F2” = (Fmeck / Poh) ee, where F, obs is the mock 
observational flux, and F™°% is the mock flux of a galaxy 
SED in each of the five WFST bands. 

Dust extinction is taken into account when generating mock 
flux data. The SED flux density after dust reddening from 
interstellar medium (Calzetti et al. 1994; Galametz et al. 2017) 


can be expressed as 


Sextinct (Arest) = Sintrinsic (Arest) Jj: 9580-198 A 


(2) 


where E(B — V) = Ay/Ry is the color excess and K(A) is the 
dust extinction curve. We adopt the Calzetti et al. (2000) 
extinction curve, with Ry for this attenuation law set as 4.05. 
For each galaxy, the value of attenuation Ay is given by the 
COSMOS/UltraVISTA catalog, which is derived through the 
SED fitting technique. We directly use it to generate the mock 
extinction-corrected fluxes. 

We also consider intergalactic medium (IGM) absorption for 
high-redshift galaxies. At wavelengths shorter than the Lya 
line, the emission can be absorbed by neutral hydrogen clouds 
in the IGM along our line of sight to the high-redshift galaxy. 
We account for this extinction by making use of the Madau 
(1995) IGM attenuation law. This is carried out by applying the 
average flux decrement (D4) between Lya and Ly, and (Dy) 
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Figure 3. Some best-fit galaxy SED examples in the observed frame in the COSMOS/UltraVISTA field (Muzzin et al. 2013). 


between Lyf and the Lyman limit, such that the IGM 
absorption corrected flux can be written as 


Sabsorption (Arest) — e Y (Di)) Sinitial (Arest) (i = A, B), (3) 


where Sjnitiai is the initial flux density in the rest frame, adopted 
as the interstellar dust extinction-corrected galaxy flux 
Sextinct(Arest) Obtained from Equation (2). After these correction 
procedures, the galaxy SED flux density Sybsorption, With dust 
extinction and IGM absorption corrected, is substituted into 
Equation (1) to generate mock flux data for all five WFST 
bands. 

Next, we estimate flux errors with respect to mock WFST 
fluxes. For a ground-based telescope, the signal to noise ratio 
(S/N) can be evaluated via the following equation (Lei et al. 
2023), 


S/N 
= S-A-+T 
SS-A- T +2 > ny. (Sky A- og + D) T +R] 
(4) 


where S is the source signal with a constant spectral flux, 7 is 
the exposure time, A is the effective area of the WFST primary 
mirror (~4.12 x 10* cm?), apix = 0.111 arcsec? is the area 
of one pixel, D is the dark current of the CCD 
(D —0.005e^ pixel! s ', at — 100 ?C), R? is the readout noise 
of the CCD (R= 8e rms), and ni, is the total pixel number in 
the PSF region. The factor of 2 applied here is because we 
assume that the calculation is performed on sky subtracted 
images. We adopt an optimal PSF aperture of 1.18 times the 


full width at half maximum for a non-adaptive optics case 
according to the Integration Time Calculator of Gemini. Sky in 
Equation (4) is the sky background signal that actually lands on 
the detector in units of e^ s! pixel !, which is given by 


Sky — J FT Tani OEGsdA, (5) 


where f, is the surface brightness of the sky background, Topt is 
the throughput of the optics (including the primary mirror, 
analog to digital converters and the five corrector lenses), and 
QEccp is the quantum efficiency of the CCD. 

The photometric error can be evaluated through the 
magnitude error given by the approximate relation 
Opn © 2.5log(1 + 1/S/N) (Pozzetti et al. 1998; Bolzonella 
et al. 2000). We also add a systematic error o,y, = 0.02 mag 
(Cao et al. 2018), and the total magnitude error is then given by 


On = Ton T E Thus we can obtain the flux error op of 
each band from o,, via error propagation. Finally, a random 
error drawn from the Gaussian probability distribution function 
(with c = gp) is added to the mock flux in each band as the 
final mock photometry. 

After computing and correcting for these mock fluxes, the 
mock observational targets obtained in the WFST shallow 
mode and deep mode are generated. In this paper, we adopt 3c 
detections to include sources into various samples as in Muzzin 
et al. (2013), i.e., galaxies with fluxes that meet the 3c depth 
thresholds of the five WFST bands (see Table 1) are selected as 
the mock observational samples for subsequent Zpnot calculation 
(see Section 4). 
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Figure 4. Magnitude distributions of the full, training, and validation mock galaxy samples (indicated by different colors) in the (a) shallow mode and (b) deep mode, 
respectively, with the lunar phase being fixed to 90? (i.e., half moon). The color plus symbols show the medians of the distributions and the horizontal error bars 
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3. Computation of Photometric Redshifts 


In this paper, we compute Zphot Of galaxies using the mock 
WFST data and the ZEBRA photometric-redshift code 
(Feldmann et al. 2006) with default parameters unless stated 
otherwise. The main advantage of ZEBRA is that it can 
generate a new set of templates adaptive to observed galaxy 
SEDs to minimize the mismatch between observed SEDs and 
available templates. This is performed by creating a training set 
of galaxies to optimize the shape of spectral templates that can 
better match predicted galaxy colors with observed ones. We 
adopt the same set of nine initial galaxy templates (see 
Figure 2) as in Muzzin et al. (2013) for Zphot calculation using 
ZEBRA. Since we have removed all the point sources that are 
likely bright stars or AGNs in the COSMOS/UItraVISTA 
catalog, we do not include any AGN templates during our 
template fitting. 

First, we run ZEBRA in the photometry-check mode to 
identify and correct systematic errors in the photometry based 
on the maximum-likelihood algorithm. ZEBRA derives a 
simple photometric offset in each band that minimizes the 
residuals between the mock observed fluxes and that of the 
best-fit templates, with the redshifts set as the input ones (i.e., 
Zspec Or high-quality Zphot, if Zspec are not available, from Muzzin 
et al. 2013). These corrections are then applied to the mock 
WEST photometry data, and ZEBRA iterates this procedure for 
5 times to ensure that the median offset in each band converges. 

Second, we run ZEBRA in the non-template-improvement 
mode based on this photometric systematic offset-corrected 
mock catalog, using the nine initial galaxy templates shown in 
Figure 2. ZEBRA iteratively performs five logarithmic 
interpolations in the magnitude space between any adjacent 
pair of the nine templates, to generate 5 x 8 =40 templates 
added to the nine initial templates, resulting in a total of 49 
templates. 


Third, we run ZEBRA in the template-improvement mode, 
where ZEBRA transforms the discrete template space into a 
linearly continuous space, using a Karhunen—Loéve expansion 
to iteratively correct the eigenbases of a lower dimensional 
subspace through a x? minimization scheme. As a result, 
adaptive spectral templates are generated to better match the 
galaxy SEDs of the training set than the set of 49 templates. 

For each galaxy sample considered, we randomly divide all 
its galaxies into two equal halves: one half serves as the 
training set to generate a new set of adaptive templates, and 
Zphot COmputation is performed on the other half as the 
validation set based on these new adaptive templates plus the 
above 49 templates as a blind test of Zphot quality. We compare 
multiwavelength photometry and input redshifts of both the 
training and validation sets in Figures 4 and 5, respectively. We 
find that they have almost identical photometric and redshift 
properties, such that the templates generated based on the 
randomly selected training set of galaxies can be adaptive to the 
full galaxy sample, and that Zphot computation on the validation 
set as a blind test can be representative of the result for the full 
sample. 

In the template-improvement mode, ZEBRA iterates twice, 
i.e., over the redshift of 0—3 as one single bin and in smaller 
redshift bins of 0.5, to train the 49 templates based on a chosen 
training set Narrowing down the redshift bin (e.g. to 
Az = 0.2) only increases the total number of adaptive templates 
generated, but has little effect on the final Zphot results. 
Therefore, we use a total of 49 x 6 + 49 = 343 final templates 
and run ZEBRA to compute Zphot for each selected galaxy 
sample. 


4. Results and Discussion 


In this section, we show Zphor results with mock WFST data 
in the shallow and deep modes given various lunar phases (see 
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Figure 5. Distributions of input redshifts from Muzzin et al. (2013) of the full, training, and validation mock galaxy samples (indicated by different colors) in the (a) 
shallow mode and (b) deep mode, respectively, with the lunar phase being fixed to 90° (i.e., half moon). The color plus symbols show the medians of the distributions 


and the horizontal error bars indicate the 1c ranges. 


Sections 4.1 and 4.2, respectively), compare our WFST Zphot 
results with that from some recent works (see Section 4.3), and 
assess the improvement of Zphot quality with the addition of 
other data (see Section 4.4). 


4.1. Zpnor Results in the Shallow Mode 


The Zphot results with mock WFST data in the shallow mode 
are shown in Figure 6, whose left and right panels are for the 
non-template-improvement and template-improvement modes 
under various lunar phases, respectively. To evaluate Zphot 
quality, we adopt some commonly used quantities: (1) 
normalized median absolute deviation (e.g., Brammer et al. 


2008), i.e., anman = 1.48 x ( fe mein 
“input 


Zoutput — <input» with Zoutput and Zinput being the output Zphot and 
input redshifts from the COSMOS/UItraVISTA catalog 
(Muzzin et al. 2013), respectively; (2) outlier fraction foutiier, 
with outliers being defined as sources with 
[Az|/(1 + Zinu) > 0.15; and (3) bias, ie. median of 
Az/(14- Zinput) With outliers being removed. 

According to Figure 6, under various lunar phases, we have 
bias = —0.001-0.006, oNmap = 9.015-0.031, and fouttier = 
3.23%-5.19% in the non-template-improvement mode, and 
have bias = 0.000-0.006, øonmap = 0.011-0.029, and 
Soutlier = 3.7296—5.279?6 in the template-improvement mode, 
respectively. The template-improvement mode delivers smaller 
biases and oymap than the non-template-improvement mode, 
which is expected; however, the former mode provides 
comparable or even slightly larger foutier than the latter mode, 
due to misidentification of Lyman break as Balmer break or 
vice versa that is caused by the relatively limited photometry 
(i.e., only ugriz bands) although the significantly enlarged 
template set can cover the full parameter space of the observed 
galaxy SEDs. 


), where Az= 


As shown in Figure 7, Zphot quality shows some variation 
with lunar phases: Zphot quality improves as the lunar phase 
increases, with the best Zphot result achieved under the lunar 
phase of 180 deg (full moon). Two factors can influence zpno: 
quality of the selected sample under different lunar phases. One 
is the lunar phase itself: under brighter lunar phases, the sky 
light background contributed by the moon becomes larger, 
resulting in larger uncertainties on photometry and eventually 
Worse Zphot quality. The other is the sample selection effect: 
under brighter lunar phases, only brighter sources can be well 
observed, which usually have higher-S/N photometry that 
leads to higher-quality Zphot- 

To make a more sensible evaluation of the lunar phase 
influence and consider the above two factors separately, we 
restrict the sample observable under full moon and measure 
phot under different lunar phases, with the results shown as the 
red dashed lines in Figure 7. It is clear that the lunar phase has a 
very limited influence on Zphot results of a fixed sample. 
Therefore, the variation in Zphot quality under different lunar 
phases in the shallow mode is primarily driven by the sample- 
selection effect. 


4.2. Zpnor Results in the Deep Mode 


The zygo results with mock WFST data in the deep mode 
under various lunar phases, with ZEBRA run in the template- 
improvement mode, are shown in Figure 8. Apparently, the 
inclusion of large amounts of faint galaxies significantly 
reduces Zphot quality: anman grows from 0.041 to 0.064 with 
the dimming of lunar phases; foutlier increases to 26.6% when 
there is no moon; bias —0.005, being almost constant and 
comparable to the situations of faint lunar phases in the 
shallow mode. 

In the deep mode, Zphot quality shows a stronger variation 
with lunar phases than in the shallow mode, as shown in 
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Figure 6. zn: results in the shallow mode with ZEBRA run in the non-template-improvement mode (left panels) and template-improvement mode (right panels) 
under various lunar phases (0 deg: no moon; 45 deg: 1/4 moon; 90 deg: half moon; 135 deg: 3/4 moon; and 180 deg: full moon), respectively. In each panel, blue 
dashed lines depict the boundary of Zphot Outliers, and the number of sources considered, bias, foutiier, and anman are annotated. 
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Figure 9. However, this does not mean that dimming of 
moonlight will cause Zphot quality to decrease for a fixed 
galaxy sample. When we consider the fixed sample of 
galaxies observable under full moon, we find that dimming of 
sky background caused by moonlight slightly reduces 
photometric uncertainties and thus improves Zphot quality, 
e.g., foutlier decreasing from 17.1% (full moon) to <14% (no 
moon) (see the red dashed lines in Figure 9). Thus, the 
downgrade of Zphor quality under fainter lunar phases in the 
deep mode is a direct result of the sample-selection effect, 
same as in the shallow mode. The dimmer moonlight in the 
deep mode enables the detection of fainter populations of 
galaxies, which often exhibit poorer photometry qualities; 
consequently, this leads to a continuous decrease in the 
accuracy and reliability of Zphot estimation. Therefore, we 
conclude that lunar phase only has negligible or very slight 
effects on Zphot quality for a given sample of galaxies; 
however, it can have a strong influence on sample selection, 
resulting in "apparent" variation of Zppot quality across 
different samples. 


4.3. Comparison with Other Zpnor Results 


We compare our WFST zy, results with some relevant 
works; for simplicity, we fix the lunar phase in the WFST mock 
data to 90deg (half moon) here. Figure 10 shows 
Az/( + Zinput) as a function of r-band magnitude in the 
shallow mode and deep mode (see the black contours), 
respectively. Overall, the bright sources in the shallow mode 
have much better Zpnot than the faint sources in the deep mode, 
with the scatter of Az/(1 + Zinput) Of the latter being ~2-3 
times larger than that of the former. The red curves in Figure 10 
show the average cumulative rms deviation between Zphot and 
Zspec aS à function of r-band magnitude in the SDSS survey 
early data release (using the ugriz-band photometry; Csabai 
et al. 2003), where Zphot were derived with a hybrid technique 
(empirical and template fitting methods) to calibrate galaxy 
SED templates to improve zio; quality, utilizing a training set 
of galaxies with secure Zspe We find that, at m, « 22, our 
Az/(1 + Zinput) Scatter is generally comparable to or smaller 
than that of Csabai et al. (2003). This is partly because their 
training set of galaxies are restricted to the bright population, 
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Figure 8. Similar to Figure 6, but for Zpnot results in the deep mode with ZEBRA run in the template-improvement mode under various lunar phases. 


which makes it difficult to constrain Zphot scatter toward the 
faint end. Recently, Yang & Shen (2023) estimated zp, of 
galaxies and quasars in the Southern Hemisphere DES wide 
survey based on a Bayesian analysis algorithm in the multi- 
color space, using the grizY-band photometry. We show the 
standard deviation of Az/(1 + Zinpu) Of their galaxies in the 
blue bars in Figure 10, which is comparable to our result in the 
shallow mode. 


10 


Figure 11 shows Az/(1 + Zinput) as a function of Zinput in the 
shallow mode and deep mode, respectively, in comparison with 
several other works. In general, our Az/(1 + Zinput) Shows a 
smooth distribution in each smaller redshift bin; the biases and 
scatters of our Zphot are smaller than many quoted results from 
other works up to z ~ 3. This may be because the training sets 
we use to improve the galaxy SED templates are randomly 
selected, thereby having good coverage of various galaxy 
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properties and being representative of the full galaxy sample 
(see Figures 4 and 5). However, in real observations, Zspec of 
the training sets would be mostly limited to bright sources and 
low redshifts, being difficult to well cover the full properties of 
the selected galaxy sample; in addition, the observed galaxy 
SEDs can be very different from the galaxy templates adopted 
here; therefore, a nonnegligible effect on actual biases and 
scatters of our Zphot in real observations would be expected. 

Figure 12 shows onmap and fouuie; as a function of Zinput in 
the shallow mode and deep mode, respectively, in comparison 
with the aforementioned works. Again, our Zphot results are 
overall in line with those in the literature. At z « 1.5, our Zphot 
quality is comparable to those of most other works, but not 
better than those works based on machine deep learning, e.g., 
using random forest algorithms or convolutional neural 
networks. At z > 1.5, both our Zphot results and the quoted 
results deteriorate; our foutlier is larger than the results based on 
the 5-band HSC photometry that includes the near-infrared Y 
band conducive to Zphot improvement at high redshifts. In 
contrast, our ONmap remains largely constant and acceptably 
small both in the shallow mode and deep mode, and within the 
full redshift range of 0—3 explored here. 

It is clear that at low redshifts (z < 1.5), to a certain degree, 
the machine deep learning procedures can effectively further 
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improve Zphot results compared to the traditional template- 
fiting techniques, which is usually done by applying a large 
training sample with secure Zspec and high-quality observed 
SEDs. At higher redshifts (z 2 1.5), however, such a training 
sample would become very incomplete, which makes it 
difficult to cover the full parameter space of all observed 
sources; thus, it is still unlikely to precisely constrain 
uncertainties of Zphot measurement at high redshifts simply 
based on machine learning. At z > 1.5, the traditional template- 
fitting technique still shows advantages in some respects, e.g., 
as shown in Figure 12, our oymap outperform that of machine- 
learning results, because ZEBRA can extend the known 
templates in the multi-parameter space and improve the fitting 
result by creating new templates and optimizing their shapes to 
be adaptive to galaxy multiwavelength photometry. However, 
the ZEBRA template-improvement procedure does not seem to 
effectively reduce foutier at z 2 1.5, mainly due to misidenti- 
fication of spectral breaks or other spectral features in galaxy 
SEDs thanks to the limited ugriz-band photometry. In contrast, 
the most recent machine learning methods based on the Direct 
Empirical Photometric (DEmP) or Nearest Neighbor (NNPz) 
method seem to have the potential to reduce foutier to a large 
extent. Therefore, in the future, we can combine the machine 
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deviation of DES galaxies based on grizY-band photometry (Yang & Shen 2023). 
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Figure 11. Az/(1 + zinpu) as a function of Zinput in the shallow mode (left) and deep mode (right), respectively. Also shown for comparison are those from the HSC 
survey (using convolutional neural network for Zphot computation, NetZ; cyan; Schuldt et al. 2021), second Red-Sequence Cluster Survey (using Direct Empirical 
Photometric method, DEmP; green; Hsieh & Yee 2014), HSC-SSP survey (using K Nearest Neighbor, KNN; blue; Zou et al. 2022), DES survey (using KNN; brown; 
Zou et al. 2022), DESI survey (using KNN; purple; Zou et al. 2022), HSC survey (using DEmP; orange; Tanaka et al. 2018; no error bars provided), HSC survey 
(using Nearest Neighbor, NNPz; gray; Tanaka et al. 2018; no error bars provided), and SDSS survey (using random forest regression; pink; Carliles et al. 2010). 


learning methods with adaptive template fitting procedures to 
further improve WFST Zpnot quality. 


4.4. Improvement of Zpnor Quality with the Addition of 
Other Data 


We further investigate the improvement of WFST Zphot 
quality by including mock data from the China Space Station 
Telescope (CSST, to be launched around 2024; Zhan 2011) 
and Euclid space telescope (launched in 2023 July; Laureijs 
et al. 2012), both of which can provide additional high-quality 
ultraviolet and/or near-infrared data in large sky areas that are 
critically supplementary to WFST data. 
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We consider the CSST NUV- and y-band mock data, whose 
photometric errors are measured via S/N (Ubeda 2011): 


Cyt 


» (©) 
(Gt En Nix (Bay zx Baet)t ae Nox Nea Rz 


S/N — 


where t is the exposure time and N,jx is the number of detector 
pixels covered by a source. Npix is 16 by default, corresponding 
to the case of a point source in the image; changing N,;, value 
does not significantly alter the final result. N,eaq is the number 
of detector readouts, Bye is the detector dark current, and R, is 
the read noise. Default parameter settings of t= 300s, 
Naa = 2, Baa = 0.02e-s-^!pixel!, and R,—5e pixel! 
are adopted. C, is the count rate from the source in units 
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of e sl. Bay in Equation (6) is the sky background in 
es ! pixel '. For more details about the CSST mock flux and 
error estimation, we refer readers to Section 2.3 in Cao et al. 
(2018). 

We consider the Euclid Yg-, Jg- and Hg-band mock data. 
Since we do not have specific details of Euclid (such as those 
of CSST shown in Equation (6)), we adopt photometric errors 
in the similar Y, J, and H bands of the VISTA survey for 
approximation, i.e., photometric errors of mock Euclid data 
are directly taken from the Muzzin et al. (2013) catalog, 
which are scaled proportionally to mock Yg-, Jg- and Hg-band 
fluxes. Given that there is a slight bias between the ground- 
based VISTA telescope and Euclid, we apply a constant 
conversion factor to convert the VISTA errors to the Euclid 
mock errors, which is defined as the ratio of flux error 
between the CSST y band and VISTA Y band for each source 
at the given magnitude. We then compute mock fluxes and 
flux errors in the CSST NUV, y and Euclid Yg, Jg and Hg 
bands, which are subsequently combined with WFST mock 
data for Zphot improvement. 
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Figure 13 shows the Zphot results in the deep mode with the 
addition of 5-band mock data from CSST and Euclid. It is clear 
that the zyyo, quality is significantly improved (see Figure 8), 
because the 10-band mock photometry that well covers the 
wavelength from ultraviolet to near infrared is vital for both 
ZEBRA photometry-check mode and template-improvement 
mode. In the non-template-improvement mode, foutlier and 
Onmap ate effectively reduced to ~5% and ~0.03, respec- 
tively; the lunar phase has little influence on Zphot results, 
mainly due to the fact that mock CSST and Euclid data are 
almost unaffected by the lunar phase. In the template- 
improvement mode, foutier and ONMaAD are further reduced to 
~1% and ~0.02, respectively; meanwhile, the bias is also 
better calibrated, being ~0.0. 

Fulfillment of many scientific goals is heavily dependent on 
Zphot accuracy, €.g., Zphot for future photometric weak lensing 
surveys need to at least achieve oymap < 0.05, with many 
relevant studies setting oNmap~0.02 as a goal (e.g. 
Zhan 2006; LSST-Collaboration et al. 2009), which is crucial 
to depict the redshift dependent growth of dark matter 
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Figure 13. Similar to Figure 6, but for Zphot results in the deep mode with the addition of mock data from the CSST-NUV, CSST-y, Euclid- Yg, Euclid- Jg, and Euclid- 


Hg bands. 


fluctuations, analyze weak lensing cosmic shears, and inves- 
tigate the redshift dependent weak lensing signal behind 
clusters of galaxies under the constraints of the framework of 
the dark energy state equation (Brimioulle et al. 2008). As 
shown above, such requirements on Zphot accuracy can be met 
when the mock WFST, CSST and Euclid are combined. 


5. Summary 


In this paper, we conduct a preliminary study that assesses 
Zphot Quality based on the mock WFST ugriz-band photometry 
in the shallow mode and deep mode. We adopt the multi- 
wavelength photometric catalog in the COSMOS /UltraVISTA 
field to generate mock WFST data, as it has deeper limiting 
magnitudes than WFST observations; during this process, 
mock fluxes are computed through the convolution of galaxy 
SEDs with the five WFST filter transmission curves, with 
interstellar dust extinction and IGM absorption taken into 
account, and mock flux errors are evaluated through the 
consideration of instrumental parameters, sky background, and 
systematic errors. 
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We calculate Zphot using the ZEBRA code, which can 
generate new adaptive templates that better describe observed 
galaxy SEDs. We find bias <0.006, ewwAp S 0.03, and 
fouttier [5% in the shallow mode and bias 720.005, 
anman © 0.06, and foutlier 7 17%-27% in the deep mode, 
respectively, under various lunar phases; a lunar phase has 
limited influence on Zphor results, and the decrease of Zphot 
quality with dimming of the lunar phase is primarily caused by 
sample-selection effect, i.e., the involvement of increasingly 
more fainter sources that have larger photometric uncertainties. 

We compare our WFST Zpnor results with that of some 
relevant works, finding general agreement between various 
results. Given that the adaptive template fitting and machine 
learning methods have their respective merits, it would be 
sensible to use all these methods jointly to further improve 
WEST Zphot quality in the future. 

Finally, we compute Zphot by combining the mock WFST 
data with ultraviolet and near-infrared data from CSST and 
Euclid. As expected, we find significant improvements in Zphot 
quality with fouaie; ~ 1% and anman 7 0.02, thanks to the full 
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wavelength coverage from ultraviolet to near-infrared. Such 
high-quality Zphor can help fulfill many scientific goals that 
highly rely on Zphot accuracy. 
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